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Abstract 

We study thermal leptogenesis in a broad class of supersymmetric S'O(IO) models with a left- 
right symmetric seesaw mechanism, taking into account flavour effects and the contribution 
of the next-to-lightest right-handed neutrino supermultiplet. Assuming Mc = M u and a 
normal hierarchy of light neutrino masses, we show that four out of the eight right-handed 
neutrino mass spectra reconstructed from low-energy neutrino data can lead to successful 
leptogenesis with a reheating temperature in the (10 9 — 10 10 ) GeV range. In the remaining 
four solutions, leptogenesis is dominated by N2 decays, as in the type I seesaw case. We 
find that some of these spectra can generate the observed baryon asymmetry for reheating 
temperatures above 10 10 GeV, in contrast to the type I case. Together with flavour effects, 
an accurate description of charged fermion masses turns out to be a crucial ingredient in 
the analysis. 
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1 Introduction 



Leptogenesis is one of the most popular mechanism for generating the observed [T] baryon asymmetry 
of the universe: 

n B~ n B = (6,21 ±0,16) x 10- 10 . (1) 

In its simplest version [2], out-of-equilibrium decays of heavy Majorana neutrinos generate a lepton 
asymmetry which is then partially converted into a baryon asymmetry by sphaleron processes [31 0] . This 
mechanism has been extensively studied in the last decade [5] . In particular, conditions for a successful 
leptogenesis have been obtained [HI El E] and many refinements have been added, such as spectator 
processes [9] , finite temperature corrections [8] and flavour effects [10] . One of these conditions is the 
famous Davidson-Ibarra bound [6] on the lightest right-handed neutrino mass, M\ > 0(1O 8 — 10 9 ) GeV, 
which applies in the case of a hierarchical right-handed neutrino mass spectrum. The main outcome of 
these studies is that thermal leptogenesis can work with parameters consistent with the type I seesaw |11] 
interpretation of neutrino oscillation data. By contrast, standard electroweak baryogenesis fl2j fails to 
produce the observed baryon asymmetry |131 114] . and its supersymmetric version |15] is successful only 
in a small portion of the Minimal Supersymmetric Standard Model (MSSM) parameter space^J 

While thermal leptogenesis can successfully generate the baryon asymmetry of the universe if we 
are free to choose the right-handed neutrino masses and couplings (modulo the constraints coming from 
neutrino masses and mixing), this might not be the case if the seesaw mechanism is embedded into 
a more fundamental theory - typically a Grand Unified Theory (GUT) based on the £0(10) gauge 
group [T7J. In such theories, the right-handed neutrino parameters are constrained both by the unified 
gauge symmetry, which implies relations among quark and lepton mass matrices, and by neutrino 
oscillation data, with no guarantee that they fall into the range preferred by leptogenesis. It is well 
known indeed [18] that the 50 (10) mass formula = M u leads to a strongly hierarchical heavy 
neutrino mass spectrum (except for special values of the light neutrino parameters pjj] ) , with M\ lying 
below the Davidson-Ibarra bound. This conclusion can be evaded, however, if the relation Md = M u 
receives large corrections from Yukawa couplings involving a 126 or a 120 Higgs representation, as in the 
so-called minimal ,50(10) model [20J and its extensions, or from non-renormalizable interactions |21j . 
The contribution of the next-to-lightest right-handed neutrino to the baryon asymmetry could also, in 
principle, change the above picture [22J, |23] . 

In this paper, we investigate another possibility to reconcile 5*0(10) unification with leptogenesis, 
based on the left-right symmetric seesaw mechanism^] In a broad class of 50(10) models, neutrino 
masses receive contributions from both the type I [11] (right-handed neutrino exchange) and the type 
II [27] (heavy scalar SU{2)l triplet exchange) seesaw mechanisms, with both contributions related by 
a left- right symmetry. As a result, for a given Dirac mass matrix, 8 different right-handed neutrino 
mass spectra are consistent with the same light neutrino mass matrix [28^ , instead of a single one in the 
type I case. It was shown in Ref. [29 (albeit with a qualitative discussion of the washout) that some of 
these spectra can lead to successful leptogenesis even if the mass relation Md = M u holds. This was 
confirmed, for the case of an inverted light neutrino mass hierarchy, in Ref. [30] , where the Boltzmann 
equations were solved in the one-flavour approximation (the possibility of triplet leptogenesis has also 
been studied in Ref. [31] ) . The purpose of the present paper is to perform a more comprehensive study 
of thermal leptogenesis in this class of 5O(10) models, including previously missing ingredients such 
as flavour effects and the contribution of the next-to-lightest right-handed neutrino, and investigating 
the dependence of the final baryon asymmetry on low- and high-energy parameters as well as on the 
reheating temperature. The necessary corrections to the GUT-scale mass relation = M e are also 

2 Electroweak baryogenesis could however still be a viable mechanism in other extensions of the Standard Model in 
which the dynamics of the electroweak phase transition is modified (see e.g. Ref. [TB] for a review). 

3 A completely different option, based on a non-standard embedding of the Standard Model matter fields into 5*0(10) 
representations, has been explored in Refs. [241 125] . Successful leptogenesis can also be achieved by adding chiral singlets 
to supersymmetric SO(10) models [26] . 
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taken into account in our analysis. Assuming Mp = M u and a normal hierarchy of light neutrino 
masses, we find that successful leptogenesis is possible with a reheating temperature in the (10 9 — 10 10 ) 
GeV range for 4 out of the 8 reconstructed right-handed neutrino spectra. Some of the remaining 4 
spectra can also generate the observed baryon asymmetry, but for higher reheating temperatures. 

The paper is organized as follows. In Section [2] we review the left-right symmetric seesaw mechanism 
in supersymmetric SO(10) models, as well as the properties of the associated right-handed neutrino 
mass spectra. In Section [3j we write the flavour-dependent Boltzmann equations governing leptogene- 
sis, including the contribution of the next-to-lightest right-handed neutrino and of its supersymmetric 
partner. In Section [4] we solve numerically the Boltzmann equations and present our results for the 
final baryon asymmetry, taking into account the corrections to the mass relation = M e . In Sec- 
tion |5j we study the dependence of the final baryon asymmetry on the yet unmeasured light neutrino 
parameters, on the high-energy Dirac couplings and on the reheating temperature. Finally, we present 
our conclusions in Section HJ 



2 The framework 

2.1 The left-right symmetric seesaw mechanism in supersymmetric 50(10) models 

One of the appealing features of S'O(IO) unification is that it provides a natural realization of the 
seesaw mechanism, thus yielding an elegant explanation for the smallness of neutrino masses. Indeed, 
the right-handed neutrinos needed for the (type I) seesaw mechanism belong to the 16j representations 
(i = 1,2,3) that contain the Standard Model matter fields, and they acquire heavy Majorana masses at 
the scale where the B — L symmetry (which is part of the 50(10) gauge symmetry) is broken. Namely, 
the Majorana mass matrix is generated either by the renormalizable operators 16jl6jl26 or by the 
non-renormalizable operators 16jl6jl6 16/A, in which the S77(5)-singlet component of the 126 (resp. 
16) Higgs representation provides the (B — L)-breaking vev. As for the Dirac neutrino mass matrix, it 
arises from the same SO (10) operators that contribute to the charged fermion masses. Upon integrating 
out the heavy Majorana neutrinos, one obtains the well-known type I seesaw mass formula: 

AfW = -MlM^M D , (2) 

where and Mr denote the Dirac and Majorana mass matrices, respectively. In supersymmetric 
SO{\SS) models with a 126 Higgs representation, the light neutrino mass matrix can receive an additional 
type II contribution if a 54 Higgs representation is also present. The role of the 54 is to induce a coupling 
between the 126 representation, in which the SU(2)l triplet needed for the type II seesaw mechanism 
lies, and a 10 representation containing a significant H u component where H u is the MSSM Higgs 
doublet responsible for up quark masses. To see how this works, it is convenient to use a left-right 
symmetric language: the 126 contains a right-handed triplet A c with quantum numbers (1,1,3)_2 
under SU(3)c x SU(2) l x SU(2)r x U(l) b-l, whose vev vr is responsible for the breaking of the B — L 
symmetry, as well as a left-handed triplet A = (1,3, l)+2; the 54 contains a bitriplet A = (l,3,3)o; 
and the 10 contains a bidoublet = (1,2, 2)o- The superpotential terms relevant for the type II seesaw 
mechanism read: 

l - LiLjA + l -a $$A + r AA C A , (3) 



where the first term comes from the 16;16jl26 couplings, while the second and third terms come from 
the 1010 54 and 54126126 couplings, respectively. Setting ((A c )°) = vr in Eq. ^ and integrating 
out the heavy triplets A and A, one obtains: 

M {II) = *u}% f (4) 

2M A J ' v ' 



4 In general, H u is a linear combination of all Y = +1 SU(2)l doublets contained in the Higgs representations of the 
model. Denoting by the Y = +l doublet lying in the 10 under consideration, one can write = a u H u + • • • , where 
the dots stand for heavy Y — +1 doublets. We assume here that a u ~ 1, although strictly speaking only a u ^ is required. 
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where a u = a^cr, v u = (H®) = vsinft (v = 174 GeV), and Ma is an effective SU(2)l triplet mass. The 
couplings of Eq. Q alone would give Ma = tvr, but the superpotential generally contains additional 
terms contributing to the SU(2)l triplet mass matrix. Depending on these, Ma may be larger or 
smaller than vr (for vr <C Mgut, & tuning of the superpotential parameters is generally necessary to 
achieve Ma < vr [32]). Notice that, due to the left-right symmetry embedded in the SO (10) gauge 
symmetry, the same set of parameters fij determine the triplet couplings in Eq. Q and the right- 
handed neutrino mass matrix, which is given by Mr = /vr. Assuming further that the Dirac neutrino 
mass matrix is symmetric (which excludes a contribution from Yukawa couplings involving a 120 Higgs 
representation), one ends up with the left-right symmetric seesaw mass formula: 

M v = ^ff-^Yvf-'Y,,, (5) 
2 Ma Vr 

where we have written the Dirac mass matrix as Md = Y u v u . 

Definite predictions for the baryon asymmetry generated via leptogenesis require the knowledge 
of the masses and couplings of the heavy decaying states. In this respect, 50(10) models provide 
a predictive framework, since the Dirac mass matrix is generated by the same Yukawa couplings as 
the charged fermion mass matrices. In the type I seesaw case, one can reconstruct the right-handed 
neutrino mass matrix Mr from the knowledge of the light neutrino mass matrix by simply inverting 
Eq. ([2]), provided that the Dirac mass matrix is known. Assuming e.g. that Md and M u only receive 
contribution from renormalizable Yukawa couplings to 10- dimensional Higgs multiplets, which implies 
the well-known mass relation Md = M u , one generically obtains a strongly hierarchical right-handed 
neutrino mass spectrum, with M\ lying below the Davidson-Ibarra bound [6] (see however Ref. |19] 
for special situations where M\ and M2 can be degenerate). More generally, successful leptogenesis is 
difficult to achieve in 50(10) models with a type I seesaw mechanism, even taking into account the 
contribution of the next-to-lightest right-handed neutrino [29 as suggested in Ref. [23] . 

In the left-right symmetric seesaw case, the reconstruction of the matrix / that determines both 
the right-handed neutrino mass matrix and the triplet couplings requires the resolution of the non- 
linear matrix equation (|5]). In Ref. [28], Akhmedov and Frigerio showed that this equation has exactly 
2 n solutions in the n generation case, and provided explicit solutions up to n = 3. An alternative 
reconstruction procedure, which employs complex orthogonal matrices, was proposed in Ref. [29]. There 
it was argued, based on a qualitative discussion of the washout, that this multiplicity of solutions 
makes it possible for leptogenesis to be successful in 50(10) models with a left-right symmetric seesaw 
mechanism. The purpose of the present paper is to put this statement on a quantitative basis, and to 
prove in particular that successful leptogenesis is indeed possible for "mixed" solutions in which neither 
the type I not the type II seesaw contribution dominates in the light neutrino mass matrix. 



2.2 Reconstruction procedure 

Before presenting our study, let us briefly recall the reconstruction procedure of Ref. [29] . Our starting 
point is the left-right symmetric seesaw formula ([5]), in which both / and Y v are complex symmetric 
matrices. We want to reconstruct / for a given pattern of light neutrino masses and lepton mixing, 
assuming that Y u is known in a basis in which the charged lepton mass matrix is diagonal. For con- 
creteness, we work in the 3-family case, but the procedure applies to any number of neutrino families. 
In order to solve Eq. (|5]), we first rewrite it as 

Z = aX-(3X-\ (6) 

with a = a u vl/(2M A ), = v\jv R and 

Z = N-'M^N-y, X = N-'fiN-y, (7) 
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where N v is a matrix such that Y v = N U N^ , and Y v is assumed to be invertible. Being complex and 
symmetric, Z can be diagonalized by a complex orthogonal matrix if its eigenvalues (i.e. the roots of 
the characteristic polynomial det(Z — zl) = 0) are all distinct: 

Z = O z Biag{z 1 ,z 2 ,z 3 )0^ , \zi\ < \z 2 \ < \z 3 \ , O z O T z = 1. (8) 

ThenEq. @ can be solved for X in a straightforward manner, by noting that X is diagonalized by the 
same complex orthogonal matrix as Z. Upon an Oz transformation, Eq. ^ reduces to 3 independent 
quadratic equations for the eigenvalues of X: 

Zi = axi - fix' 1 . (9) 

For a given choice of (x±, x 2 , x 3 ), the solution of Eq. §5§ is given by: 

f = N u O z Diag{x 1 ,x 2 ,x 3 )O^N^ . (10) 

The right-handed neutrino masses Mi = j%VR are obtained by diagonalizing / with a unitary matrix: 

/ = U f Diag (fx, f 2 , h) Uj , h<h<h, U f Uj = 1 , (11) 

where the fi are chosen to be real and positive. The matrix Uf relates the original basis for right-handed 
neutrinos, in which Y v is symmetric, to their mass eigenstate basis. It can be used to express the Dirac 
couplings in terms of charged lepton and right-handed neutrino mass eigenstates, as A = U^Y V . 

Since each equation z% = ax% — (3x~ l has two solutions x~ and xf, there are 2 3 = 8 different solutions 
for the matrix /, which we label in the following way: (+,+,+) refers to the solution ( 
(+, +, — ) to the solution (xf, x 2 , £3 ), and so on. It is convenient to define x^ and xf such that, in the 
4a(3 <C \zi\ 2 limit: 

x~ ~ -j , xf ~ ^ . (12) 

We will refer to x^ as the "type I branch" and to xf as the "type II branch". This terminology is 
motivated by the fact that solutions (—,—,—) and (+,+,+) reduce to the "pure" type I and type II 
cases in the large vr limit (defined by 4a/3 <C \zi\ 2 , or equivalently vr ^> 2a u v^/ Ma\zi\ 2 ): 

/<"•-•-> 4 ^ 112 - Vl Y u M~^, (13) 

VR 

f(+,+,+) 4a/3<,!| 2 Mu_ _ 

a 

The remaining 6 solutions correspond to mixed cases where the light neutrino mass matrix receives 
significant contributions from both types of seesaw mechanisms. In the opposite, small vr limit ( | ^3 1 2 -C 
4a/3), one has xf ~ ± sign(Re(zj)) \f~j3ja, which indicates a partial cancellation between the type I 
and type II contributions to light neutrino masses. Finally, in the region of intermediate vr values, 
|zi| 2 < 4a/3 < I ^3 1 2 , both the type I and the type II seesaw mechanisms give significant contributions 
to the light neutrino mass matrix. 



2.3 Properties of the reconstructed right-handed neutrino spectra 

The above procedure can be used to determine the a priori unknown couplings in theories which 
predict the Dirac matrix Y u , taking low-energy neutrino data as an input. In Ref. &$\, it was applied 
to supersymmetric .50(10) models with two 10's, a 54 and a pair of 126 © 126 representations in the 
Higgs sector, but no 120 representation (as required by the left- right symmetric seesaw mechanism). In 
this subsection, we recall the main properties of the reconstructed right-handed neutrino mass spectra. 
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Given the above assumptions, the most general renormalizable Yukawa couplings read: 

16*16^10! + yA 2) 16*16^102 + fij 16*16,126 , (15) 
where Y^ and / are complex symmetric matrices. Assuming that the SU(2)l doublet components 



of the 126 do not acquire a vev, Eq. (15) leads to the following GUT-scale mass relations: 



M u = M D , M d = M e . (16) 

It is well known that the second relation is in conflict with experimental data and needs to be corrected. 
In the absence of a 210 Higgs representation that would induce vev's for the doublet components of 



the 126 [33; , this must be done by non-renormalizable interactions. We postpone the discussion of this 
issue to Section [4] and assume for the time being that Eq. (16) holds. 



The inputs in the reconstruction procedure are the matrices Y u and M u at the seesaw scale. Ne- 



glecting the running of Y v between the GUT scale and the seesaw scale, Eq. (16) yields, in the basis for 



the 16 matter representations in which M e is diagonal with real positive entries: 

Y v = UjY u U q , U q = PuVcKMPd, % = Biag (y u ,y c ,y t ) , (17) 
where Vqkm is the CKM matrix and y u c ,t are the U P quark Yukawa couplings, all renormalized at 



the GUT scale. The presence of two diagonal matrices of phases P u and Pd in Eq. (17) is due to the 
fact that the £0(10) symmetry prevents independent rephasing of right-handed and left-handed quark 
fields. In the same basis, the light neutrino mass matrix generated from the seesaw mechanism reads: 

M v = UfMJj} , Ui = PeVpMNsPu , M v = Diag (mi, m 2) m 3 ) , (18) 

where Vpmns contains a single, Dirac-type phase Spmns, and P e and P v are two diagonal matrices 
of phases. With the convention that P v contains only two phases, Vpmns = VpmnsPv is the PMNS 
lepton mixing matrix and mi^,3 are the light neutrino masses, all renormalized at the seesaw scak£J The 
two phases in P v are the physical CP-violating phases associated with the Majorana nature of the light 
neutrinos, while the three phases contained in P e , analogous to the five independent phases contained 
in P u and Pd, are pure high-energy phases. Once the input values for (y u ,yc,yt), (mi,m2,m3), Vckm 
and Upmns at the GUT scale are fixed, the 8 different / matrices can be reconstructed as a function of 
a, (or, equivalently, of the B — L breaking scale vr and (3 /a) and of the high-energy phases contained 



in P u , Pd and P e . Notice that, as long as Eq. (16) holds, the reconstructed fij couplings depend on 



the combination of phases PdP e rather than on Pd and P e separately; hence the number of independent 
high-energy phases reduces to five. 

Let us specify the input values that we are going to use in this paper. For the quark masses and the 
CKM parameters at the M z scale, we take the central values given in Refs. [M] and [35], respectively: 

m u {M z ) = 1.7 MeV, m c (M z ) = 0.62 GeV, m t {M z ) = 171 GeV, (19) 
m d {M z ) = 3.0 MeV, m s {M z ) = 54 MeV, m b {M z ) = 2.87 GeV, (20) 
A{M Z ) = 0.818, A = 0.2272, p = 0.221 , a- 0.340. (21) 

These values, together with the inputs for the lepton mass and mixing parameters, are subsequently 
evolved to the GUT scale M GUT ~ 2 x 10 16 GeV using the Mathematica package REAP [36] with 
an effective supersymmetric threshold MgusY = 1 TeV and tan/3 = 10. Our reference light neutrino 
spectrum is a normal hierarchical spectrum with m\ = 10 -3 eV and #13 = 0. For the observed oscillation 
parameters, we take the best fit values of Ref. |37j : 

Am| 2 = m\-m\ = 2.4 x 10~ 3 eV 2 , sin 2 9 23 = 0.44, (22) 
Amjjj = m\-m\ = 7.92 x 10~ 5 eV 2 , sin 2 12 = 0.314. (23) 



5 Strictly speaking, Eq. (|5J) involves the decoupling of four states at scales which can differ by several orders of magnitude. 
We neglect the associated radiative corrections here and, for simplicity, identify the seesaw scale with the GUT scale in 
the following. 
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The Dirac-type phase Spmns an d the two relative Majorana phases contained in P v are treated as 
free parameters, as well as the high-energy phases contained in P u , P^ and P e . We denote these 
phases by q>^ d,u ' e : % = 1,2,3 (some of which are redundant), with P u = Diag (e 4 *i , e**2 , e l *3), P d = 
Diag (e**i , e**2 ; e **3) 5 and so on. 

Fig.[T]shows four out of the eight right-handed neutrino spectra reconstructed from the above inputs 
as a function of vr, assuming (3/a = 0.1. The 8 different solutions can be distinguished by the behaviour 
of each Mi (i = 1,2,3) as a function of vr. The 4 solutions with 23 = x~^ are characterized by a constant 
value of the lightest right-handed neutrino mass, Mi ~7x 10 GeV; among them, the 2 solutions with 
x 2 = %2 a ^ so have M2 ~ 4 x 10 9 GeV, while the 2 solutions with 22 = x t have a rising Mi. It is 
interesting to note that 7 x 10 4 GeV and 4 x 10 9 GeV are nothing but the type I values of M\ and M2, 



respectively. This is actually not a coincidence, but a direct consequence of Eq. (12): besides the fact 
that solution (—,—,—) reduces to the type I case in the large vr limit, some of the properties of the type 
I right-handed neutrino mass spectrum are inherited by the four solutions^] (db, db', — ). The 2 solutions 
with X3 = and X2 = x^, on the other hand, are characterized by M\ ~ 2 x 10 9 GeV, and the 2 
solutions with £3 = x% and X2 = x\ by a rising M\. In the large vr region, solution (+, +, +) approaches 
the type II case, as can be seen from the fact that M\ : M2 : M3 (oc /1 : /2 : fz) oc mi : rri2 : m^. It 
should be stressed that the choice of (3/a does not affect the shape of the curves Mj = M{(vr), but 
only their position along the horizontal axis. For instance, setting (3/a = 1 instead of 0.1 would shift 
the curves in Fig. [I] according to vr — > VO.l vr. This is due to the fact that, while / depends on a and 
(3 separately, Mr = Jvr only depends on the combination a(3. In the following we choose (3/a = 0.1 
for numerical convenience, but it should be clear that our results will not depend on this choice. 

The implications for leptogenesis of these solutions were discussed in Ref. |29j. Since Ma = 
((3/a) u u vr/2 < ((3/a)vR and all solutions satisfy Mi <C vr, one can safely assume that the SU(2)i 
triplet is heavier than the lightest right-handed neutrino. Then the dominant contribution to the lepton 
asymmetry comes from the out-of-equilibrium decays of N\ and Ni, except in some cases where the 
contribution of N2 and A^2 actually dominates. One can distinguish between three different behaviours 
(from now on, Ni will refer both to the ith right-handed neutrino and to its supersymmetric partner): 

• solutions (d=, dz', — ): the four solutions characterized by a low value of Mi fail to generate the 
observed baryon asymmetry from JVi decays. In these solutions, the CP asymmetry in N\ decays 
always lies below 0(1O~ 10 ), irrespective of the choice of the high-energy phases. This situation is 
similar to the one encountered in the type I case. In principle, A^2 decays could generate a large 
asymmetry in a lepton flavour that is only mildly washed out by N\ decays and inverse decays, 
thus leading to successful leptogenesis [23] • Estimates of this effect tend to show, however, that 
it is unlikely to work in solution (—,—,—). 

• solutions (=b, +, +): in the two solutions characterized by a rising Mi, the CP asymmetry in N\ 
decays grows with vr. Successful leptogenesis then becomes possible for large values of vr, i.e. 
in the region where the type II seesaw contribution dominates in the light neutrino mass matrix. 
However, one finds a tension between successful leptogenesis, which requires Mi > 10 10 GeV, and 
the gravitino overproduction problem |38| . which imposes an upper bound Trjj < (10 9 — 10 10 ) 
GeV on the reheating temperature [39] . 



solutions (dz, — ,+): the two solutions characterized by Mi ~ 10 9 GeV can lead to a relatively 
large CP asymmetry in N\ decays without conflicting with the gravitino constraint. However, 
the washout of the generated lepton asymmetry by lepton number violating processes tends to be 
large. To determine whether the observed baryon asymmetry can indeed be generated, one must 
integrate numerically the Boltzmann equations. 



6 A comment about our notation might be necessary: the prime symbol in ±' means that the second sign in (±, ±', — ) 
is not correlated with the first one, i.e. (±, ±', — ) refers to the four solutions (+, +, — ), (+, — , — ), (— , +, — ) and ( — , — , — ). 
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Case + ,+ ,+ 



Case + ,-,+ 



1C) 16 
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Figure 1: Right-handed neutrino masses as a function of vr in solutions (+,+,+), (+,—,+), (+,—,—) and 
( — ,—,—). Inputs: hierarchical light neutrino masses with mi = 1CP 3 cV, oscillation parameters as specified 
in the text, and no CP violation besides the CKM phase (6 PM ns = = $f = = $? = 0); f3/a = 0.1. 
The range of variation of vr is restricted from above by the requirement that \fij\ < 1. Dashed lines indicate a 
cancellation at a stronger level than 1% between the type I and type II contributions to the light neutrino mass 
matrix. 



A general feature of all solutions is that lepton number violating processes tend to efficiently wash out 
the generated lepton asymmetry. This can be traced back to the relation Mjj = M u , which implies that 
at least one of the Dirac couplings is of the order of the top quark Yukawa coupling. As a consequence, 
predictions for leptogenesis depend on the details of the dynamics encoded in the Boltzmann equations. 

It is clear from the above discussion that, for most solutions, the qualitative analysis of Ref. [29] 
is not sufficient to tell whether leptogenesis can indeed be successful. The purpose of the present 
paper is to perform a careful, quantitative study of leptogenesis in supersymmetric SO (10) models with 
a left-right symmetric seesaw mechanism, taking into account the lepton flavour dynamics [40 -|53|. 
as well as the contribution of the next-to-lightest right-handed neutrino supermultiplet |22| [23] [52} 
|54"1 155] . As is well known from studies performed in the type I case, flavour effects can significantly 
enhance the final baryon asymmetry if there is a hierarchy between the washout parameters for different 
lepton flavours [HI H2J H3], and their impact might be crucial for solutions (±, ±', — ) and (±, — ,+). 
Furthermore, the contribution of the next-to-lightest right-handed (s)neutrino can be relevant both for 
Mi ~ M.2 and for M\ <C M2, provided in the latter case that iVi-related washout effects are weak. The 
Boltzmann equations including all these ingredients will be presented in the next section. Finally, the 
corrections to the GUT-scale mass relation = M e needed to account for the measured down quark 
and charged lepton masses will modify the reconstructed seesaw parameters and affect the final baryon 
asymmetry. These effects will be taken into account in Section [4] 
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3 Boltzmann equations 



In this section, we write the Boltzmann equations that govern thermal leptogenesis in the class of 
supersymmetric SO (10) models described above. The relevant heavy degrees of freedom are the three 
right-handed neutrinos iVi ; 2,3 and the scalar Higgs triplet A, as well as their supersymmetric partners. 
In principle, all these states could contribute to the generation of the lepton asymmetry through their 
decays. However, due to the strong hierarchy among their masses, only the decays of N\ and N2 and 
of their scalar partners are relevant in practice^} 

Let us first consider the CP asymmetries in heavy (s)neutrino decays. The relevant superpotential 
terms, written in the basis of charged lepton and right-handed neutrino mass eigenstates, read: 

W seesaw = \ a NfL T a icj 2 H u + l -MiNfNZ + ^f a/3 Lli<T 2 AL + ~a u H^ia 2 AH n + M A Tr(AA) , (24) 

where N?, L a and H u are the right-handed neutrino, lepton doublet and Higgs doublet superfields, 
respectively, and 

Being Majorana particles, the heavy neutrinos can decay both into £ a H u (lepton + Higgs boson), £ a H u 
(slepton + higgsino) and into the CP-conjugated final states £ a H* and £* a H u . Their scalar partners 
Nf are not CP eigenstates and have only 2 two-body decay modes, Nf — ► £ a H u and Nf — ► £ a H u . The 
corresponding tree-level decay rates are [56 : 

M- 

T(N t - £ a H u ) = T(Ni - £ a H*) = T(N { - £ Q H U ) = T(N{ - t a H u ) = — - \X ia \ 2 , (26) 

167T 

T(N9 -> £ a H u ) = T(Nf -> £ a H u ) = ^ |A iQ | 2 . (27) 
Supersymmetry implies the equality of the total decay widths, = T-^- c = M^XX^h/Att. At the 

i 

one-loop level, an asymmetry between CP-conjugated decay channels arises from the interference of 
tree- level and one- loop Feynman diagrams. One can define 4 different types of flavour-dependent CP 
asymmetries which, due to supersymmetry, are given by the same quantities ej Q : 



tic* — 



T(Ni - £ a H u ) - T(Ni - £ a H*) T(N t - £ a H u ) - T(N, - t a H u ) 



Y(N i - £H U ) + T(N l - £H*) T{Ni ^ m u ) + T(Ni - £*H U ) 

r(N9 - 4# u ) - r(ivf - r(ivf - l a H u ) - v(Wf - # *) 



T(N?^£H U )+T(N? -^£H U ) T(Nf ^ £H U ) +T(Nf £*H*) 



(28) 



where r(JV< -» lff„) = E^iV; -» etc. The €iaS receive both a type I contribution (right- 

handed neutrino-induced vertex and self-energy corrections) and a type II contribution (triplet-induced 
vertex correction): 

£ia = + (-ia ■ (29) 

The type I contribution reads [56 : 

1 Im[A tQ (AA%A*J (Mf\ 

7 Indeed, all solutions satisfy Mi,2 <C M3,Ma over a large range of values for vr (assuming Ma ~ ([3/a)vn and /3/a 
not too small), hence the lepton asymmetry generated in A and ^3 decays is washed out by JVjj- and iVi-related processes. 
In the large vr region of some solutions, one has instead M2 ~ M3 (Ma), but N2 and N3 (A) turn out to be too heavy 
to be thermally produced after reheating. Indeed, one typically has M2 > 10 12 GeV in this case, while the upper bound 
on the reheating temperature associated with the gravitino overproduction problem is O(10 9 — 10 10 ) GeV. 



8 



while the type II contribution is given by (571 [58] : 



ji 



8tt 



E 



Im[Ai /3 (Mj / )* Q A iQ ] Mi ( M 



(AAt)j 



/// 



m 2 



(31) 



In Eqs. (30) and (31), the loop functions are given by: 

2 



1 — x 



In 



1 + x 



(32) 



(33) 



Mi ~ M2 <C M3, they must be modified. Following Ref. |59| . we take: 



Eqs. (30) and (32) are valid for hierarchical right-handed neutrinos; in the case of a partial degeneracy, 

(34) 



1 1 

^(aaT),, ^ 

3^1 



^ Im (X ia X* ja 



(G? + 2C^) + (AA^Cg,}) , 



where Cj 7 , C|f a and C* J . are functions of M|/M?. is the usual vertex function: 



x In 



while C^ a and C/. arise from self-energy corrections: 



1 + x 



X 



(35) 



x C^, (x 



x(l — x) 



(1 - x) 2 + x (AAt)l /167T 



(36) 



The term proportional to C l J b in Eq. (34~| does not contribute to the total CP asymmetry in iVj decays 
and is therefore a pure flavour effect [5TJ . In the limit of hierarchical right-handed neutrinos, x> 1, 
this term can be neglected and C^ a (x) ~ yfx/(l — x). One then recovers the non-resonant formula (30). 

The strength of the processes involving Ni {Nf) and £ a (£ a ) is parametrized by individual washout 
parameters^] Ki a : 



T(Ni 



y H u ) + T(Ni 



Defining the effective neutrino masses 



H(Mi[ 



\Xia\ 2 vl 

Mi 



the individual washout parameters be expressed as 

^ia — j 



(37) 



(38) 



(39) 



where m* = lGn^^y/g^v^/ (3\/5Mp) ~ (1.56 x 10 _3 eV) sin 2 j3 is the equilibrium neutrino mass (we used 
the fact that the number of effectively massless degrees of freedom in the thermal bath at T = M\ is 
g* = 228.75 in the MSSM). Summing over flavour indices, one obtains the total washout parameters Kf 



(1.56 x lO" 3 eV) sin 2 (3 



(40) 



We follow here the notations and conventions of Refs. 
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The out-of-equilibrium condition for JVj decays (resp. JVj decays into £ a or i a ) reads «j < 1 (resp. 
Ki a < 1). If JVj is in the strong washout regime (/tj S> 1), scattering processes in the thermal bath 
produce an equilibrium population of JVj at T ~ Mj. When T drops below Mj, the JVj decay and 
generate asymmetries in all lepton flavours, in proportions controlled by the ej Q; a = e, r. Assuming 
that the decays occur in the temperature regime where all flavours are relevant, the dynamical evolution 
of the individual flavour asymmetries depends on the Kj a , where the index j runs over the JVj such that 
Mj < Mj. If Kj a S> 1, the asymmetry in the lepton flavour a is strongly washed out by the lepton 
number violating processes involving Nj , which are in equilibrium at T = Mj ; if Kj a <C 1 , the opposite is 
true. Consider now the case where JVj is in the weak washout regime (Kj <C 1). Then its number density 
never reaches thermal equilibrium, but since K{ a <C 1 for all a, the individual flavour asymmetries 
are only weakly washed out by JVj-related processes (they may however be strongly washed out by 
JVj-related processes, j / i). 

The evolution of the number densities is obtained by solving the set of Boltzmann equations. As is 
usual done, we take into account the expansion of the universe by defining comoving number densities 
Yx = nx/s. The supersymmetric Boltzmann equations for the heavy (s)neutrino comoving number 
densities read 1611: 



Y^{z) = -2 Kl (D l (z) + S l (z)){Y Ni (z)-Y N! ' ! i (z)) , (41) 
Y' n {z) = -2« i (AW + 5 i (z))(y^ 4 (z)-l^(z)) , (42) 
where Y^. stands for Y^ c + Y~j , z = M\/T and the symbol ' stands for d/dz. The equilibrium densities 



appearing in Eqs. (41) and M2h are given by: 



yjfW = ti?z 2 K 2 (Riz) ^ !M aUx 10 -3 , (43) 

where Ri = Mi /Mi, and we have corrected the high temperature behaviour of the Maxwell-Boltzmann 
distribution by a factor of 3£(3)/4 in Y^ 9 (z), and by a factor of £(3) in Y^(z) |^J The factor of 



2 in the right-hand side of Eqs. (41) and (42) accounts for the fact that there are twice as many 
channels in the supersymmetric case as in the non-supersymmetric case. In the regime in which all 
lepton flavours are relevant, the individual A a = B/3 — L a asymmetries are driven by the following 
Boltzmann equations 



Y' Aa {z) = -2Y,^^m^) + S l {z))[Y Ni {z)-Y^{z) + [Y^{z)-Yg{z))) 

i=l,2 

+ 2 Ki^Wiiz) A aP Y A p(z) , (45) 

i=l,2 

where Y Aa stands for the total A a asymmetry stored in the fermionic species and in their supersym- 



metric partners. In Eqs. (41), (42) and (45), the thermally averaged decay rates Di(z) are given by: 



9 Assuming Maxwell-Boltzmann (MB) statistics as is customary, the equilibrium number density at high temperature 
differs from the Fermi-Dirac (FD) and Bose-Einstein (BE) cases by a numerical factor: 

n FD (T » M) = \ n BE (T » M) = ^ n MB (T » M) , 
while n FD (T) ~ n BE (T) ~ n MB (T) at low temperature (T < M). 



10 



The scatterings terms Si{z) account for Higgs-mediated AL = 1 scatterings involving top quarks and 
antiquarks. They receive both s- and t-channel contributions: 



Si(z) = 2Sl(z)+4Si(z) , (47) 

whose expression can be found in Ref. [60]. The washout term Wi(z) = W^ D (z) + Wf(z) results from 
the contribution of inverse decays: 

w! D (z) = ^Rf z 3 Ki(Riz) , (48) 

and AL = 1 scatterings [60 : 



^m-^(^>(^ + ^j + "*m>- (49) 

In writing the above Boltzmann equations, we made several assumptions which we now proceed 
to clarify. In the washout term, we neglected the off-shell part of the AL = 2 scatterings, which is a 
good approximation as long as Mj <C Ki a (10 13 GeV) |43j. We also omitted AL = scatterings such as 
NiNj — ► ££, N{Nj — > H U H* and N{£(£) — > Nj£(£), which do not contribute to the washout but can affect 
the abundance of the heavy (s)neutrinos (when flavour effects are taken into account, they also tend to 
redistribute the lepton asymmetry among flavours). These processes are of higher order in the neutrino 
Yukawa couplings and are expected to have little impact on the final baryon asymmetry. We further 
neglected the triplet-related washout processes, gauge scatterings [8] [59], spectator processes |9j, and the 
higher order processes 1 — > 3 and 2^3 [62]. Finally, since the left-right symmetry is broken at a scale 
vr which may lie several orders of magnitude below Mqut, one may worry that decay and scattering 
processes mediated by the heavy SU(2)r x U(1)b-l gauge bosons affect the heavy (s)neutrino number 
densities. Indeed, Wr- and Z'-mediated processes such as Ni en — > Ijr q' R and NiNi — ► // tend to keep 
the heavy (s)neutrinos in thermal equilibrium, thus reducing the generated lepton asymmetry |63|, 164] . 
As shown in Ref. [BS], however, this effect can be practically neglected if M^/vr < 10~ 2 , which turns 
out to be the case for i = 1, 2 in each of the 8 solutions (at least as long as Ni % are light enough to be 
thermally produced after reheating, i.e. M\ t 2 % 10 11 GeV). Therefore, we do not need to include Wr- 
and Z'-mediated processes in our stud} 
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The Boltzmann equations ( 45 ) are written for the L?/3 — L a asymmetries 1a q rather than for the 



lepton asymmetries Y^ a , because the former are preserved by all MSSM interactions (including the 



non-perturbative sphaleron processes), contrary to the latter. As the washout term in Eq. (45) depends 
on Yf, a , we need to express it in terms of the Ya^s. This is done by a conversion matrix A [TU], whose 
entries depend on which interactions are in equilibrium: 

Y La = J2 A ^ Y ^ ■ ( 5 °) 

& 

Depending on the temperature at which leptogenesis takes place (identified for simplicity with M\ 
below), one must consider one of the following three regimes. If M\ < 10 9 GeV (1 + tan 2 /?), the tau and 
muon Yukawa interactions are in equilibrium, hence all three lepton flavours are distinguishable. One 
must then write three separate Boltzmann equations for the flavour asymmetries Y& e , Ya„ and Yj\ T , as 



in Eq. ( 45 ) . The 3 x 3 A matrix is given by [45 : 



-93/110 6/55 6/55 
A = | 3/40 -19/30 1/30 | . (51) 
3/40 1/30 -19/30 



10 Note however that Wr - and Z'-mediated scatterings can help generating an equilibrium population of Ni in the weak 
washout regime (k,: 1), provided that Mi/vn > 10 -3 |65j . 
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If 10 9 GeV(l + tan 2 /?) < Mi < 10 12 GeV(l + tan 2 /3), the muon Yukawa interactions are no longer 
in equilibrium, and the electron and muon lepton flavours can no longer be distinguished. The lepton 
asymmetry must then be projected onto the 2-flavour space (Yl b+ , Yl t ) , where ^ / L e+fi = Yl b + Y^^, 
and correspondingly Y/± e+ii = Y/± e + Y"a m • The conversion matrix is now a 2 x 2 matrix: 

/ -541/761 152/761 \ 
" V 46/761 -494/761 J ' 1 > 

and the Boltzmann equations for Ya b and Ya m must be replaced by a single equation for Ya e+m involving 
the CP asymmetries ej je + M = £j e + and the washout parameters Hi, e +n = Kie + k^. Finally, if 
Mi > 10 12 GeV (1 + tan 2 none of the interactions involving charged lepton Yukawa couplings are in 
equilibrium, and one recovers the flavour-independent treatment of leptogenesis with A = — 1. 

It has been pointed out in Ref. (HJ] that the out-of-equilibrium condition used above to determine 
the flavour regime is actually not sufficient. For a given flavour asymmetry Y& a to evolve independently 
during leptogenesis, the corresponding charged lepton Yukawa interaction rate T a (T) ~ 5 x 1CP 3 /i 2 (1 + 
tan 2 (3) T should not only be faster than the expansion rate of the Universe, but, more importantly, it 
should also be faster than the Ni inverse decay rate. If this were not the case, inverse decays would 
keep the evolution of the lepton state coherent. For the tau and muon lepton flavours, this condition 
reads Mi < 10 12 GeV (1 + tan 2 [3)/k 1t and Mi < 10 9 GeV (1 + tan 2 (3) / m^, respectively. It is more 
restrictive than the out-of-equilibrium condition for K\ T 3> 1 (resp. fti„ 3> 1). 

The final baryon asymmetry is given by: 



1 D 



« 



where the factor 10/31 is due to the partial conversion of the A Q asymmetries into a baryon asymmetry 
by the non-perturbative sphaleron processes [1] (we assume here that sphalerons come out of thermal 
equilibrium below the electroweak phase transition). It has been pointed out recently [66 that the 
conversion factors relating the flavour asymmetries Y& to Yb depend on the actual superpartner 



spectrum. Strictly speaking, Eq. (53) is only valid if all sfermions are heavy, while the final baryon 
asymmetry can be reduced by a factor of 2/3 if sleptons are light, as in some minimal supergravity 
scenarios. Since we do not assume any particular superpartner spectrum in this paper, we shall stick 



to Eq. (53) in the following. 

Being supersymmetric, the Grand Unified models we are considering in this paper face the so-called 
gravitino problem [35]: in an inflationary universe, gravitinos are abundantly produced by thermal 
scatterings in the reheating phase and, if they are unstable, their late decays tend to spoil the successful 
predictions of Big Bang nucleosynthesis. The requirement that this does not happen puts an upper 
bound on the reheating temperature, Trh ^ (10 5 — 10 10 ) GeV for 300 GeV < m 3 / 2 < 30 TeV, depending 
on the gravitino mass and on the superparticle spectrum [67] . If instead the gravitino is the lightest 
supersymmetric particle (LSP), the requirement that its relic density does not exceed the dark matter 
abundance leads to a weaker constraint, Trh < (10 9 — 10 10 ) GeV for a gravitino of mass m 5 /2 ~ 100 
GeV [39 . These constraints are in conflict with successful thermal leptogenesis, which requires Mi > 10 9 
GeV in the type I seesaw case [6]. In the left-right symmetric seesaw case studied in this paper, we also 
find such a tension. We nevertheless stick to thermal leptogenesis, since some supersymmetric scenarios 
can accommodate a reheating temperature Trh ~ 10 10 GeV (see the discussion at the end of Section[5]). 
In our numerical computation of the baryon asymmetry, we do not explicitly include the dynamics of 
reheating in the Boltzmann equations, as was done in Refs. [8] |4"H1 [68] . but we take the reheating 
temperature into account in the initial conditions. Namely, we start evolving the Boltzmann equations 
at some temperature Tj n , which we identify with Trh- Heavy (s)neutrinos with masses Mj > (4 — 5) Tj n 
will thus give a negligible contribution to the lepton asymmetry, because their production processes 
(inverse decays and AL = 1 scatterings) are Boltzmann-suppressed. 
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4 Numerical results 



In this section, we present our numerical results for leptogenesis in the class of supersymmetric SO (10) 
models described in Section |2j The final baryon asymmetry is obtained by numerically integrating 
the Boltzmann equations written in Section |3j starting the evolution at Tj n = 10 11 GeV with vanishing 
abundances for iVi 2 and iVi 2. With this choice of Tj n and tan (3 = 10, one can consider that leptogenesis 
always takes place in the 3-flavour regime. In Section [5] we shall investigate the effect of lowering T*i n . 

We present results for the following four reference solutions: 

• (+, +, +): this solution, which corresponds to dominance of the type II seesaw contribution in the 
large vr region, is characterized by a mild hierarchy of right-handed neutrino masses, with M\% 
growing with vr. 

• (+, — , +): this solution is characterized by an intermediate value of M\ {Mi ~ 10 9-10 GeV). 

• (—,—,—): this solution, which reduces to the type I seesaw case in the limit vr — * 00, is 
characterised by a hierarchical right-handed neutrino mass spectrum with a small value of Mi 
(Mi ~ 10 4-5 GeV) and an intermediate value of M2. 

• (+,+, — ): this solution differs from the previous one by the fact that M2 grows with vr rather 
than assuming a constant value. 

The four remaining solutions show very similar patterns of right-handed neutrino masses (the main 
differences are in the behaviour of M2 and M3 in the large vr region) and give analogous results for 
the baryon asymmetry. The input values for the quark masses and mixing angles and for the measured 
neutrino oscillation parameters are chosen as specified in Section [2j We further assume a hierarchical 
light neutrino mass spectrum with mi = 10~ 3 eV and #13 = 0; in Section [5] we shall study the effect 
of varying these input parameters. Rather than attempting to perform a scan over the large number 
of available CP-violating phases, we choose specific values for these phases in order to illustrate the 
typical behaviour of the different solutions. We first show the results obtained under the assumption 
that Md = M e holds at the GUT scale; in a second stage, we include the necessary corrections to 
this mass relation and investigate their effects on the right-handed neutrino mass spectrum and on the 
final baryon asymmetry. Finally, we set (3/a = 0.1 for practical reasons of numerical integration. As 



explained in Section 2.3 changing the value of (3/a only amounts to shift the curves along the vr axis. 

For convenience, we shall plot the absolute value of the baryon asymmetry \Yr\ rather than Yr 
itself. Since the sign of Yr can be reverted by simply changing the sign of all CP-violating phases, we 
do not loose any information by doing so, at least to the extent that the effect of the CKM phase (whose 
sign is determined experimentally) can be neglected. In practice this will be the case in all examples 
studied in this paper, because the CKM phase always appears in combination with small quark mixing 
angles, and large high-energy phases are assumed to be present. 

4.1 Relevance of flavour effects 

We first show in Fig. [2] the final baryon asymmetry in the absence of corrections to the GUT-scale mass 
relation M^ = M e , for the four reference solutions (+,+,+), (+,—,+), (+,+,—) and (—,—,—). In 
order to estimate the relevance of flavour effects, we plotted the result of the numerical computation 
both in the one-flavour approximation (dashed black line) and in the flavour-dependent approach (solid 
red line) . One can see that flavour effects tend to enhance the baryon asymmetry by up to one order of 
magnitude in the (+, +, +) and (+, — , +) cases. Not surprisingly, solution (+, +, +) leads to successful 



leptogenesis for large values 11 of vr, where Mi > 10 10 GeV; flavour effects allow this solution to be 



11 We do not show the region vr > 2 x 10 14 GeV, where Yb drops below the WMAP value. This behaviour is due to 
the fact that, above vr — few x 10 14 GeV, N\ and N2 are too heavy to be thermally produced (Mi : 2 3> Tin)- 
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Case +++ 



Case +-+ 




V R (GeV) V R (GeV) 

Figure 2: The final baryon asymmetry as a function of vr for the four reference solutions, in the one-flavour 
approximation (dashed black line) and with flavour effects taken into account (solid red line). The GUT-scale 
mass relation Md = M e is assumed. Inputs: hierarchical light neutrino masses with mi = eV, 6*13 = and no 
CP violation in the PMNS mixing matrix; $2 = 7r/4 and all other high-energy phases are set to zero; j3/a = 0.1. 
The Boltzmann equations are evolved starting from T in = 10 11 GeV. The thick horizontal line corresponds to the 
WMAP constraint. 



successful for smaller values of vr (i.e. for smaller values of Mi) than in the one-flavour approximation. 
By contrast, solution (+,—,+) fails to generate the observed baryon asymmetry due to the strong 
washout by inverse decays and AL = 1 scatterings, and this conclusion still holds for different choices 
of the CP-violating phases. 

Flavour effects have a much more dramatic impact in the (+,+,—) and (—,—,—) cases, which are 
characterized by a strong hierarchy between Mi and M2. In these solutions, the observed enhancement 
of Yb is due to the fact that the asymmetry in a particular lepton flavour is only mildly washed 
out by iVi-related processes, while the total washout is strong. As a consequence, the asymmetry 
generated in N2 decays is completely washed out in the one-flavour approximation, while its projection 
on this particular flavour survives when flavour effects are taken into account. This effect, which has 
been first identified in the type I seesaw framework in Ref. [23], will be discussed in greater detail in 



Subsection 4.3 Despite the huge increase in Yb, however, solution (— , — , — ) fails to reach the WMAP 
level, while solution (+,+,—) is marginally successful for vr ~ 10 14 GeV, where M2 ~ (for larger 
values of vr, M2 3> T, n and N2 no longer contributes to Yb, which then drops well below the WMAP 
value) . 
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4.2 Relevance of the corrections to the mass relation M d = M e 



The results shown in Fig. [2] were obtained assuming M d = M e to hold at the GUT scale, which only 
agrees at the order of magnitude level with the measured down quark and charged lepton masses. Let 
us now include the necessary corrections to this relation and investigate their effects on the final baryon 



asymmetry. Given the assumptions made about the Higgs sector in Section 2.3 these corrections must 



come from non-renormalizable operators 12 The simplest possibility is to add the following terms to the 
superpotential (see e.g. Ref. [69]): 

^ 16^16^-10145 , (54) 
where we assume that the Y=+l Higgs doublet in 10i does not acquire a vev, so that the mass relation 



Md = M u is left untouchec The operators (54) will modify the reconstructed couplings fij (hence 



the right-handed neutrino masses and couplings) by introducing a mismatch U m between the bases of 
left-handed charged lepton and down quark mass eigenstates: 

M\M e = Ml , M ] d M d = UlMjUm , M\M u = Ulu\MlU q U m , (55) 

where the mass matrices are written in the basis of charged lepton mass eigenstates, and M e , M d , M u 
are diagonal eigenvalue matrices. Note that M d and M e are no longer symmetric in the original basis, 



since they receive antisymmetric contributions from the operators in Eq. (54). The measured down 
quark and charged lepton masses do not determine completely the unitary matrix U m , but constrain 
its mixing angles to lie in restricted ranges (see the Appendix for details). One of them can be taken 
to be either large (#^ ~ 1) or somewhat smaller (^ ~ 0.2 — 0.3). 

Fig. [3] shows the final baryon asymmetry for four representative choices of U m . The main difference 
with the Md = M e case is that several choices for U m lead to successful leptogenesis in the solution 
(+, — ,+). This is an interesting result, since this solution is special to the left-right symmetric seesaw 
mechanism: it does not correspond to dominance of either the type I or the type II seesaw mechanism 
in the light neutrino mass matrix in the large vr limit. As for solution (—,—,—), we were unable to 
find values of U m and of the Majorana and high-energy phases allowing the final baryon asymmetry to 
reach the observed value. There is a general tendency for Yr to reach larger values for intermediate 
values of vr (typically 10 13 GeV < vr < 10 14 GeV for our choice 0/a = 0.1), where type I and type II 
contributions partially compensate for each other in the light neutrino mass matrix. 

The enhancement of Yr observed for some choices of U m can be explained by the influence of 0^ 
on the CP asymmetries ei a and on the washout parameters Ki a (see the discussion at the end of the 
Appendix). In particular, solution (+, — ,+) is found to be successful for large values of 6*^,, as in the 
sets 1 and 2 of the Appendix. We note in passing that a large mixing in U m also enhances lepton 
flavour violation, so that processes like \x — > cy or r — > ^7 might be close to their experimental limits 
(and a significant portion of the supersymmetric parameter space is already excluded). To illustrate 
the effect of U m on the right-handed neutrino mass spectrum, we also plotted Mi and M2 as a function 
of vr in Fig. [3j The enhancement of Yr with respect to the M d = M e case in solutions (+,—,+) and 
(— , — , — ) is correlated with, respectively, an increase of M\ by a factor of 10, and an increase of M2 by 
a factor of 5. One can also notice that successful leptogenesis is generally associated with right-handed 
neutrino masses above 10 10 GeV, which indicates a conflict with the gravitino problem. Some of the 
solutions that are successful for Tj n = 10 11 GeV mig ht actually fail for T in < 10 10 GeV, and we can 
already anticipate that this will be the case for solution (+, +, — ). We shall come back to this point in 
Section [5] where the dependence of Yr on the reheating temperature will be discussed. 



12 One could alternatively relax the assumptions made in Section 2.3 and introduce a 210 Higgs representation in order 
to generate vev's for the SU(2)l doublet components of the 126. In this case, the /y couplings would contribute both 
to the left-right symmetric seesaw formula and to the down quark and charged lepton masses, which would render their 
reconstruction much more difficult. 

13 If the operators 16 4 16_, 10245 were also present, they would give an antisymmetric contribution to Mu and the 
reconstruction procedure would no longer be applicable. We assume here that they are forbidden by some symmetry. 
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Figure 3: Final baryon asymmetry (left panels) and masses of Ni and N2 (right panels) as a function of vr in 
the four reference solutions with a non-trivial U m and a non-vanishing Majorana or high-energy phase. The solid 
green, dashed blue, dotted purple and dash-dotted red lines corresponds to the sets 1, 2, 3 and 4 described in the 
Appendix, respectively. The other input parameters are as in Fig. [2] 
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Tabic 1: Values of ti a and Ki a (i = 1, 2; a = e, /x, r) in solutions (— , — , — ) and (+, +, — ). The input parameters 
are chosen as in Fig. [2j and the B — L breaking scale is vr = 10 14 GeV. 

4.3 Flavour effects and N 2 leptogenesis 

Let us now discuss in more quantitative terms the interplay of flavour effects and N 2 decays in the four 
solutions (±, ±', — ), which are characterized by a strong hierarchy between the masses of the two lightest 
right-handed neutrinos. As mentioned above, the lepton asymmetry generated by N 2 is exponentially 
washed out by iVi-related processes in the one-flavour approximation. Since N\ has a small coupling 
to a particular lepton flavour, however, the asymmetry in this flavour is only mildly washed out in the 
flavour-dependent treatment, and this explains the spectacular enhancement of Yb observed in Fig. [2] 
We shall refer to this situation as "flavour-dependent N 2 leptogenesis" in the following. 

For flavour-dependent N 2 leptogenesis to be possible, some conditions on the CP asymmetries and 
washout parameters must be satisfied. In Table[T] we list the values of ti a and Ki a (i = 1, 2; a = e, //, r) 
in the (—,—,—) and (+,+,—) solutions for vr = 10 14 GeV, assuming the same input parameters 
as in Fig. [2] (in particular, U m = 1, & 2 = tt/4 and all other CP-violating phases are set to zero). 
Both solutions exhibit a similar pattern of flavoured parameters: the CP asymmetries in N± decays 
are extremelly small (ei Q < 10 -13 ), and the washout induced by N\ is strong except for the electron 
flavour (Ki e = 0.04, while K\n ~ K\ r ~ 16). By contrast, the CP asymmetries in iV2 decays are in the 
(10~ 7 - 10" 6 ) range, and the ^-induced washout is moderate. 

Let us see in detail how these features explain the results observed in Fig. [2j focusing on solution 
(+, +, — ) for definiteness (for earlier discussions of flavour-dependent N 2 leptogenesis in the type I 
seesaw framework, see Refs. [23] and [52]). N 2 decays first generate asymmetries (lA a )jv 2 in all three 
lepton flavours. Since M\ <C M 2 , the processes involving Ni are out of equilibrium at T ~ M 2 , and 
the (Y& a ) N can be computed taking into account ^-induced washout only. Neglecting off-diagonal 
entries in the A matrix, one obtains: 

O^aJ^ * -4 x 10- 10 , (Y A J N2 ~ -4 x 10- 11 , (Y At ) N2 ~ -10- 9 , (56) 

while in the one-flavour approximation the B — L asymmetry induced by N 2 is {Yb-l)n — — 6 x 10 -10 . 
As the Universe cools down, A^i-related washout processes come into equilibrium, and the evolution of 
the ^a^'s is then governed by the following Boltzmann equations: 

Y Aa (z) = 2n la A aa W 1 (z)Y Aa (z) + 2k 1oi ^ A^W^Y^z) , (57) 



in which the source term proportional to e\ a has been neglected because of its smallness. Eq. (57) yields 
the formal solution: 

Y Aa (z) ~ {Y Aa ) N2 e 2A -^^ dxW ^ +2K la Y,A aP [ Z dxW^Y^e 2 ^^^^^ , (58) 

fi+a Jzin 

where the first term corresponds to the depletion of Y Aa due to A^i-related washout processes, whereas 
the second term represents the effect of the flavour mixing induced by the off-diagonal entries in the A 
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Case (+ H — ) with off -diagonal terms. 



Case (+ H — ) without off-diagonal terms. 
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Figure 4: Evolution of the asymmetries as a function of z = M\/T in the (+, +, — ) solution, with the off-diagonal 
entries of the A matrix included (left panel) and omitted (right panel) . The thin lines represent the three lepton 
flavour asymmetries: Ya c in blue (medium grey), in green (light grey), Ya t in purple (dark grey), while the 
thick red (medium grey) and black lines stand for Yb and Y^-ti respectively. The input parameters are chosen 
as in Fig. [2j and the B — L breaking scale is vr — 10 14 GeV. 



matrix. Neglecting the off-diagonal entries of the A matrix for the moment, and omitting for simplicity 
the scattering terms in Wi(z), one obtains: 

Yi a ~ eX^»(y A J Wa , (59) 

where we have used J" °° dz z 3 Ki(z) = 3ir/2. Since K\ e C 1 C «4/j(t)> the asymmetry in the electron 
flavour is almost unaffected by iVi-induced washout, while (Ya m )jv 2 and (Ya t ) N2 are exponentially 
diluted, namely by a factor of order 10 -11 . The final baryon asymmetry is 
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Yb -^YL - ^0-92(1aX 2 * -1.2 xlO- 10 , (60) 

in good agreement with the numerical result. In the one- flavour approximation instead, the B — L 
asymmetry generated in N2 decays is completely washed out by iVi-related processes: 

YZi * e-^(YZ A L ) N2 * 6 x 10- 35 {YZi) N2 , (61) 

so that the dominant contribution to Y^ F J^ actually comes from N\ decays, in spite of the smallness of 
e\ (an analogous statement can be made about Y"a m and Ya t m the flavour-dependent treatment). All 
these results are illustrated in the right panel of Fig. |4} 

Let us now add the effect of the off-diagonal entries in the A matrix. The contribution to Ya q of the 



second term in the right-hand side of Eq. (57 ) has been evaluated in Ref. [51] , in the non-supersymmetric 
case: 

This flavour mixing does not affect Y& e , but prevents the complete depletion of Ya„ and Ya t : 

YK d (T) - 0.12 Yl e ~ -4.4 xHT 11 . (63) 

The final baryon asymmetry is only marginally affected, reaching Yg ~ —1.5 x 10~ 10 . These analytic 
estimates are confirmed by the numerical results shown in the left panel of Fig. |4j 

14 As explained earlier in this section, the sign of Yb is not relevant since it can be reverted by changing the sign of $2 
(if one neglects the small contribution of 5c km to Yb). 
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We conclude that flavour effects play a crucial role in the four solutions (±, ±',— ), where they 
render N2 leptogenesis possible. In spite of the spectacular enhancement of the final baryon asymmetry 
with respect to the one-flavour approximation, however, successful flavour-dependent N2 leptogenesis 
is difficult to achieve, as shown by Figs. |2]and[3| This is even more so in solution (— , — , — ), which in 
the large vr limit reduces to the type I seesaw case, in which flavour-dependent N2 leptogenesis was 
originally proposed as a way to achieve successful baryogenesis in GUTs [23 . 



5 Dependence on the input parameters 

The results presented in the previous section were obtained for fixed values of the light neutrino pa- 
rameters: we assumed a normal mass hierarchy with mi = 10~ 3 eV, #13 = and no CP violation in the 
PMNS mixing matrix, while the oscillation parameters were taken from Ref. [37] . Furthermore, the re- 
lation Md = M u was assumed and the Boltzmann equations were evolved, from Ti n — 10 11 GeV. In this 
section, we numerically study the influence of these input parameters on the final baryon asymmetry. 



5.1 Dependence on the light neutrino parameters 

Let us first study the impact of the yet unmeasured light neutrino parameters: the lightest neutrino 
mass, the mixing angle #13 and the Dirac phase 5pm ns, an d finally the type of mass hierarchy. 



5.1.1 Lightest neutrino mass (normal mass hierarchy) 

Since the flavour-dependent CP asymmetries e± a are bounded by (note that the upper bound is the 
same as in the type I case [13]): 



I ^ max _ ^ M\m max I Ki a 

\ei a \ < e la = ^^—J— , (64) 



and the type I inequality k\ > mi/m* does not hold, one may expect that successful leptogenesis is 
easier to achieve for quasi-degenerate light neutrinos, mi > 0.1 eV. However, varying mi also modifies 
the reconstructed right-handed neutrino parameters, which in turn affects the CP asymmetries and 
washout parameters. In particular, the right-handed neutrino masses are modified by an increase of 
mi in the following way (the right-handed neutrino couplings Aj a are also affected via the Ut matrix): 
the Mj's belonging to a "type I branch" decrease, while the ones belonging to a "type II branch" rise. 
This behaviour can be understood by noticing that, to a good approximation, the right-handed neutrino 
masses are proportional to the Xj's (see the Appendix B of Ref. |29J). Since raising mi increases the value 



of the Zi's, Eq. (12) implies that the Mj associated with some Xj ("type I branch") decreases, while the 
Mi associated with some ("type II branch") shows the opposite behavioui 15 These considerations 
explain to a large extent the impact of a variation of mi on the final baryon asymmetry, since the 
upper bound on e\ a is proportional to M\ and Ki a oc \/M\: an increase in M\ tends to enhance the 
CP asymmetry in N\ decays and to reduce the N\ washout parameters. An analogous statement can 
be made about M2 and the iV 2 -related leptogenesis parameters. 

Fig. [5] shows the region of the (mi, vr) parameter space where \Yr\ > Y^ MAP for solutions 
(+, +, +), (+, -, +) and (+, +, -), and where \Y B \ > 0.1 Y^ MAP for solution (-, -, -), which fails to 
generate the observed baryon asymmetry. The choice for U m and the high-energy phases corresponds to 
the set 1 described in the Appendix, the other input parameters being fixed as in Fig. [2} In the (+,+,+) 
case, increasing mi amounts to shift the range of vr for which leptogenesis is successful towards lower 
values. This is consistent with the fact that M\ and M2 grow with mi, and that the thermal production 



15 Strictly speaking, this is only true for 4q/3 <C \zj\ (vr S> 2a u v u /M/^\zj\ ). In the opposite limit (which is not relevant 
for the discussion below), Xj is almost independent of Zj and the associated Mi is not affected by an increase of mi. 
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Figure 5: Regions of the (mi, vr) parameter space where \Yb\ > Y^ MAP for solutions (+,+,+), (+, — ,+) and 
(+,+,—), and where \Yb\ > 0.1 Y^ MAP for solution (—,—,—). These regions are delimited by the thick black 
contour in the (+,+,+) case, the dashed red contour for (+,—,+), the long-dashed blue contour for (+,+,—), 
and the thin black contour for ( — , — , — ). Inputs: set 1 of the Appendix for U m and the high-energy phases; other 
input parameters as in Fig. [2] 



of Ni is Boltzmann- suppressed for Mi > T{ n . The behaviour of the other solutions is more interesting: 
in all three cases, the final baryon asymmetry is suppressed for a quasi-degenerate light neutrino mass 
spectrum. In the two solutions in which N2 leptogenesis plays a crucial role, namely (+, +, — ) and 
(—,—,—), this is due to the fact that the iVi-induced washout becomes strong for all flavours, as a 
result of the decrease of M\ (also, for (— , — , — ), Mi decreases with growing mi). In the (+, — , +) case, 
a larger m\ implies a smaller M\ and thus reduces the final baryon asymmetry. As a result, leptogenesis 
fails for a quasi-degenerate light neutrino mass spectrum in all reference solutions but (+, +, +). For the 
input parameters used in Fig. [5j successful leptogenesis requires mi < 0.01 eV for solution (+, +, — ), 
and mi < 0.05 eV for solution (+, — , +). 

5.1.2 0i3 and 5 PM ns 

We now turn to the dependence of the final baryon asymmetry on #13 and Spmns, the two unknown 
light neutrino parameters which control the amount of CP violation in oscillations. In Fig. [6] we first 
show the effect of varying #13 alone, from our reference value #13 = 0° to the experimental upper limit 
013 = 13°, assuming SpMNS = 0- The choice for U m and the high-energy phases corresponds to the 
set 1 of the Appendix. Furthermore, the Boltzmann equations are evolved from a somewhat lower 
initial temperature than in the previous plots: Tj n = 7 x 10 9 GeV for (+,+,+) and (+,—,+), and 
Ti n = 5 x 10 10 GeV for (+, +, — ) and (— , — , — ). The other input parameters are chosen as in Fig. [2] 
As can be seen from Fig. [6j increasing #13 generally reduces the final baryon asymmetry, especially in 
solutions (±,±,— ) where the effect is particularly pronounced. Solution (+,+,+) behaves differently, 
although in this case too the maximum value of Yb is obtained for small values of #13. 

One could be tempted to conclude from Fig. [6] that (at least for the chosen input parameters) 
successful leptogenesis favours small values of #13. However, it is not legitimate to impose Spmns = 0: 
since #13 and Spmns always appear in combination in Upmns, ° ne should study their joint effect on the 
final baryon asymmetry. This is done in Fig. [7J for the same choice of input parameters as in Fig.[6j but 
for a fixed value of the B — L breaking scale vp. One can see that successful leptogenesis is compatible 
with a "large" value of 6*13 (#13 > 5°) as soon as 5pmns is allowed to be different from zero. Such values 
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Figure 6: The final baryon asymmetry as a function of vr in the four reference solutions, for Spmns — an d 
different values of 9 13 : 6 13 = 0° (black), 2° (purple), 5° (blue), 9° (red) and 13° (green / light grey). Inputs: 
set 1 of the Appendix for U m and the high-energy phases; T in = 7x 10 9 GeV for (+, +, +) and (+, — , +), while 
Tin = 5 x 10 10 GeV for (+, +, — ) and (— , — , — ); other input parameters as in Fig. El 




Figure 7: Contour lines of the ratio \ Yg\/Y^ MAP in the four reference solutions, as a function of #13 and Spmns- 
The input parameters are the same as in Fig. [gJ and the B~L breaking scale has been fixed at vr — 5 x 10 13 GeV 
for (+,+,+) and (+,-,+), and at n fl = 6x 10 13 GeV for (+,+,-) and (-,-,-). 
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Figure 8: Region of the (m 3 , vr) parameter space where \Yg\ > Y^map m ^ e ^_|_^ _|_ } _|_^ solution, assuming an 
inverted light neutrino mass hierarchy (delimited by the black contour). Inputs: set 1 of the Appendix for U m 
and the high-energy phases; other input parameters as in Fig. [2] 



of #i3 are within the reach of upcoming reactor and first generation superbeam experiments (for a brief 
review, see e.g. Ref. [70J. For instance, #13 just below the present experimental limit together with 
SpMNS ~ 57r/8 is compatible with Yb = Y^ MAP both in the (+, +, +) and in the (+, — , +) solution. 
Although Fig. [7] was obtained for a specific set of input parameters, we can conclude that successful 
leptogenesis is possible for values of and Spmns such that CP violation in the leptonic sector can 
be established in future neutrino superbeam experiments. 



5.1.3 Inverted mass hierarchy 

Assuming an inverted light neutrino mass hierarchy leads to significantly different results from the 
normal hierarchy case, although the gross qualitative features are preserved (in particular, solution 
(+,+,+) generally leads to successful leptogenesis, while for solution (+,—,+) this depends on the 
choice of U m and of the Majorana and high-energy phases). In this paper, we do not attempt to 
perform a general study of the inverted hierarchy case, which has already been investigated in Ref. |30| . 
in the one- flavour approximation and assuming M<j = M e . To illustrate some of the differences with the 
normal hierarchy case, we just display in Fig. [8] the plot analogous to the one in Fig. [5} Only solution 
(+, +, +) is represented there, since for the choice of input parameters made in Fig. [5]none of the other 
reference solutions leads to successful leptogenesis in the inverted hierarchy case. This is already a 
noticeable difference with the normal hierarchy case. As far as solution (+,+,+) is concerned, some 
differences with respect to the normal hierarchy case can be seen in the region where the lightest neutrino 
mass lies below a few 10 -3 eV. We were not able to reproduce the result of Ref. [3D], which finds that 
solution (+, — , +) can generate the observed baryon asymmetry in the absence of corrections to the mass 
relation = M e , in the region where M\ ~ M%. This may be due to the differences in the treatment 
of leptogenesis: the study of Ref. [30J was performed in the one-flavour approximation, assuming initial 
thermal abundances for the right-handed neutrinos, while we solved the flavour-dependent Boltzmann 
equations with Ti n = 10 11 GeV and vanishing initial abundances. Furthermore, Ref. [30] used the non- 



supersymmetric analogue of Eq. (32 ) in the computation of the CP asymmetry, while we used Eq. (34), 
which reproduces the correct behaviour of the self-energy contribution in the limit of exactly degenerate 
right-handed neutrinos. 



22 



Case Case H — l- 




16 



Vfi (GeV) V R (GeV) 

Figure 9: The final baryon asymmetry as a function of vr for different values of t/2, from y2/yc(MGUT) = 0.1 
(yellow/light grey) to y2/Vc{M GUT ) = 10 (blue/dark grey). The reference case y 2 = Vc is plotted in black. Left 
panel: solution (—,—,—), set 1 for U m and the high-energy phases; right panel: solution (+,—,+), set 4 for U m 
and the high-energy phases. The other input parameters are as in Fig. [2] 



5.2 Impact of corrections to the mass relation Mr, = M u 

In the preceding subsections, we studied the dependence of the final baryon asymmetry on the values 
of the yet unmeasured light neutrino parameters. Let us now turn to the influence of the high-energy 
Dirac couplings. So far we assumed that the mass relation Md = M u holds at the GUT scale, while 
Md = M e receives corrections from non-renormalizable operators. In this subsection, we study the effect 
of a departure from Md = M u . More specifically, we assume that Md and M u are still diagonal in the 
same basisp^j but that their eigenvalues differ (yi ^ y Ui )- This has a direct impact on the right-handed 
neutrino mass spectrum, since the Mj associated with some x~ is to a good approximation proportional 
to y| in the regime vr 3> 2a u v^/M^\zj\ 2 , while the M, associated with some xt is independent of yj 
(see the Appendix B of Ref. [29 J. In particular, one has M\ oc y\ in solution (+, — , +) and M2 oc y\ 
in solution (—,—,—). One thus expects that raising 2/2 will enhance the final baryon asymmetry by 
increasing the ei a 's in the former case, and the £2 a 's in the latter case. 

This is shown in Fig. [9] in which (y2/2/c)(-^Gcr) is varied between 0.1 and 10 in solutions (+, — , +) 
(right panel) and (—,—,—) (left panel). We can see that the final baryon asymmetry increases with 
growing y2 in both solutions. In particular, successful leptogenesis becomes possible in the (—,—,—) 
case for large enough y2 (for 2/2 = Wy c , however, N2 becomes too heavy to be thermally produced 
above vr ~ 10 14 GeV, which results in the Boltzmann suppression of Yr). This conclusion is however 
dependent on the input Tj n = 10 11 GeV: it does not hold for the more realistic choice Tj n = 10 10 
GeV (see the discussion in the next subsection about the gravitino problem). In the (+,—,+) case, 
successful leptogenesis is possible for values of vr as large as a few 10 16 GeV, and this conclusion also 
holds for Ti n = 10 10 GeV. This is an interesting result, since gauge coupling unification favours a one- 
step breaking of the 50(10) symmetry, with a B — L breaking scale close to the GUT scale (a lower 
B — L breaking scale is however not excluded [7T]). Fig. [9] also shows that y2 > y c allows solution 
(+, — ,+) to be successful with a U m containing only small mixing angles (set 4), thus alleviating the 
constraints on the superpartner spectrum coming from the non-observation of lepton flavour violating 
processes such as fi — ► e'y. 

16 This is a natural assumption if the CKM matrix mainly comes from the down quark sector. In this case, and in 
the absence of cancellations between the different contributions to Mjj and M u , both matrices have a strong hierarchical 
structure with mixing angles smaller than the CKM angles. The relative rotation between the bases in which Md and M u 
are diagonal can then be neglected in the reconstruction procedure. 
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Figure 10: Regions of the (vr, T in ) parameter space where \Yb\ > Y^ MAP for solutions (+,+,+), (+, — ,+) and 
(+,+, — ), and where \Yg\ > 0.1 Y^ MAP for solution (— , — , — ). These regions are delimited by the thick black 
contour in the (+,+,+) case, the dashed red contour for (+,—,+), the long-dashed blue contour for (+,+,—), 
and the thin black contour for (— , — , — ). Inputs: set 1 of the Appendix for U m and the high-energy phases; other 
input parameters as in Fig. [2] 



5.3 Dependence on the reheating temperature 

The numerical results presented so far were obtained starting the evolution of the Boltzmann equations 
at Ti n = 10 11 GeV, in the approximation where the dynamics of reheating is neglected. In this ap- 
proach, Ti n can be identified with the reheating temperature. In order to estimate how severe the tension 
between successful leptogenesis and the gravitino problem is, we therefore proceed to study the depen- 



dence of the final baryon asymmetry on T{ n . Fig. |10| shows the regions of the {vr^ parameter space 



where \Y B \ > Y B VMAP for solutions (+,+,+), (+,-,+) and (+,+,-), and where \Y B \ > 0.1Y B VMAP 
for solution (—,—,—). The choice for U m and the high-energy phases corresponds to the set 1 of the 
Appendix, the other input parameters being fixed as in Fig. U One can see that solution (+,—,+) 
succeeds in generating the observed baryon asymmetry for values of Ti n as low as 5 x 10 9 GeV, whereas 
solutions (+,+,+) and (+,+,-) require T in > 7 x 10 9 GeV and T in > 3 x 10 10 GeV, respectively. 
While these numbers have been obtained for a particular choice of the input parameters, they unam- 
biguously show that successful leptogenesis can be achieved with a reheating temperature below 10 10 
GeV in solutions (+, +, +) and (+, — , +). As for solution (+, +, — ), Tj„ > 10 10 GeV was found to be 
a necessary condition for successful leptogenesis for all sets of input parameters we considered. This 
allows us to conclude that, for generic input parameters, the solution (+,+,—) fails to generate the 
observed baryon asymmetry if the reheating temperature is lower than 10 10 GeV. 

As discussed at the end of Section |3j there are strong constraints on the reheating temperature from 
gravitino cosmology, and this potentially conflicts with successful thermal leptogenesis. Nevertheless, 
some supersymmetric scenarios can accommodate a reheating temperature in the (10 9 — 10 10 ) GeV range, 
as required for solutions (+, +, +) and (+, — , +) to generate the correct amount of baryon asymmetry. 
One possibility is that the gravitino is the LSP; the constraint that its relic density does not exceed 
the dark matter abundance reads Trh ^ (10 9 — 10 10 ) GeV for m 3 / 2 ~ 100 GeV [39] . This scenario 
is further constrained by the requirement that the NLSP decays do not alter the success of Big Bang 
nucleosynthesis (BBN), which can be satisfied e.g. by a sneutrino NLSP [72] or by assuming some 
amount of .R-parity violation |73j. Another way of avoiding the strong constraints on the reheating 
temperature is to assume an extremely light gravitino |74j . m.3/2 < 16 eV [75] (where the upper bound 
comes from WMAP and Lyman-a forest data), or a very heavy gravitino [76 , m.3/2 > 50 TeV. In 
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the former case, the gravitino decouples when it is still relativistic and escapes the overproduction 
problem, while the NLSP is sufficiently short-lived to decay before BBN. In the latter case, the gravitino 
decays before nucleosynthesis and does not affect the light element abundances; furthermore, the LSPs 
produced in its decays are within the observed dark matter abundance for Trh < 10 10 GeV [77]. 



6 Conclusions 

In this paper, we studied thermal leptogenesis in a broad class of supersymmetric ,50(10) models with 
a left-right symmetric seesaw mechanism, including flavour effects and the contribution of the next-to- 
lightest right-handed neutrino supermultiplet. Assuming Mjj = M u and a normal hierarchy of light 
neutrino masses, we found that successful leptogenesis is possible with a reheating temperature in the 
(10 9 — 10 10 ) GeV range for 4 out of the 8 reconstructed right-handed neutrino mass spectre 17 corre- 
sponding to solutions (±, +, +) and (±, — , +). In the remaining 4 solutions, leptogenesis is dominated 
by N2 decays, as in the type I seesaw case; among those, solutions (±,+, — ) succeed in generating 
the observed baryon asymmetry for reheating temperatures above 10 10 GeV. These results show that 
50(10) models in which the Dirac mass matrix and the up quark mass matrix have a similar hierarchical 
structure are compatible with successful thermal leptogenesis if the seesaw mechanism is of the left-right 
symmetric type. As a byproduct, we found that solution (— , — , — ), which mimics the type I seesaw in 
the large vr limit, fails to generate the right amount of baryon asymmetry in this limit. This suggests 
that successful flavour-dependent N2 leptogenesis in ,50(10) models with a type I seesaw mechanism 
requires large corrections to the mass formula Md = M u , especially if one insists on Tre < 10 10 GeV. 

Both flavour effects and the corrections to the mass relation M^ = M e were found to be crucial 
ingredients for the success of solutions (±, — , +) and, for Trh > 10 10 GeV, of solutions (±, +, — ). In 
the former case, flavour effects increase the final baryon asymmetry by up to one order of magnitude, 
while in the latter case they allow a particular lepton flavour asymmetry generated in N% decays to 
survive the washout by iV"i-related processes, resulting in a spectacular enhancement of the final baryon 
asymmetry with respect to the one-flavour approximation. We also studied the dependence of the 
results on the unknown light neutrino parameters (#13, Spmns an d the type of mass hierarchy), as 
well as on corrections to the mass relation Mo = M u and on the reheating temperature. Moderate 
deviations from Md = M u were shown to extend the region of parameter space in which leptogenesis is 
successful; in particular, solution (+, — , +) was found to be successful for values of the B — L breaking 
scale as large as 10 16 GeV, as preferred by supersymmetric gauge coupling unification. 

We believe that the results presented in this paper (which extend and put on a more solid basis 
the ones of Refs. |29t 150] ) make S'O(IO) models in which neutrino masses originate from the left-right 
symmetric seesaw mechanism more attractive. Other aspects of the phenomenology of these models, 
such as lepton flavour violation (which was briefly discussed in Ref. [29 ), could discriminate further 
between the 8 different right-handed neutrino mass spectra, and are worth studying in detail. We 
note in particular that the typical values of the mixing angles in U m allowing solution (+, — , +) to be 
successful tend to enhance the branching ratios of processes like \x — > and r — > ^7. 

Finally, it would be interesting to investigate how the results presented in this paper would be modi- 
fied if the SU (2) l doublet components of the 126 were allowed to acquire a vev. In such a case, the mass 
relations Md = M e and Mo = M u receive corrections proportional to the matrix /, and the measured 
down quark and charged lepton masses can be accounted for without introducing non-renormalizable 
interactions. However, a technical complication arises from the fact that the reconstruction procedure 
cannot be performed independently of the charged fermion mass fit. We defer the study of this case to 
future work. 



17 Although we only presented results for four reference solutions, we checked that the solutions that differ by xf 
such as (+, — , +) and (— , — , +), show similar behaviour for Yb- 
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A Appendix: corrections to the mass relation Md = M e 

In this appendix, we discuss the corrections to the GUT-scale mass relation = M e arising from 
non-renormalizable operators of the form: 

^ 16,16,(10! x45) |120 , (A.l) 

where only the Y = — 1 SU{2)l doublet in 10i acquires a vev, so that the mass relation Mp = M u 
is not affected, and the 10i and 45 Higgs representations are contracted in an effective 120 represen- 
tation, implying that the couplings Kij are antisymmetric. Such operators generate a contribution to 
the down quark and charged lepton masses when the 45 acquires a vev in the T3/2 or B — L direc- 
tion. Indeed, the decomposition of the tensor product (10 x 45)| 12 o under the Pati-Salam subgroup 
SU(2)l x SU(2)r x 577(4) c contains a (2,2,1) representation generated by (2,2,l)i x (1,3,1)45, 
and a (2, 2, 15) representation generated by (2, 2, l)io x (1, 1, 15)45. The Clebsch-Gordan coefficients 
needed to distinguish from M e arise from (2, 2, 15). 

The most general situation occurs when the two vev's ((1, 3, 1)) in the T^r direction and ((1, 1, 15)) 
in the B — L direction belong to different 45's. We denote the scale of these vev's by 1*3 and vi§ and 



assume that they have GUT scale values. The other dimensionful parameter appearing in Eq. (A.l) 
is A, which we identify with the scale at which the unified gauge coupling becomes non-perturbative, 
A ~ 10 Mgut [29] • The corrected mass matrices read: 



M d = M]° + f-| «i + V ~f *») , M e = Mj° + (-| Kl - 3 V -f k 2 ) v?> , (A.2) 



where is the contribution of the 10's, v d 01 is the vev of the Y =—l Higgs doublet in 10i, and the 
matrices k\ and K2 contain the non-renormalizable couplings associated with the two 45's. 

In order to study the corrected mass matrices, it is convenient to switch to the basis where the 
symmetric contribution M}° is diagonal. Then and M e can be parametrized as: 



M d = Si fJ, 2 £3 , M e = XlSl [12 -X3E 3 , (A.3) 






with in real and Ej, X{ complex (//j and are dimensionful parameters). This structure simplifies in 
some cases. If the vev's in the T3R and B — L directions are carried by the same 45, there is a single 
matrix of non-renormalizable couplings k and all Xi are equal. In the absence of a vev in the T3R 
direction, one has x\ = X2 = ^3 = 3. The matrices M^M^ and MlM e can easily be diagonalized in the 
case of hierarchical entries, \x\ -C [i 2 "C ^3 and |ei| <C | ^2 1 *C | £"3 1 , yielding: 

m 2 b ~ ^ + 2|e 3 | 2 , (A.4) 

m 2 T ~ ^ + 2|x 3 e 3 | 2 , (A.5) 

m h m s ~ |^2/^3 +£3 1 , (A.6) 

M2/X3 + a;|e||. (A.7) 
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The fact that £3 appears in all four relations (A.4)-(A.7) implies that antisymmetric corrections can 
at most accommodate a GUT-scale r — b mass difference of the order of m^, and this conclusion also 
holds for more general hierarchies of the \Xi and £j parameters. This is not enough to account for the 
value of (m T — m&) obtained by running the measured down quark and charged lepton masses from 
Mz to Mqut = 2 x 10 16 GeV (assuming an effective supersymmetric threshold Msusy = 1 TeV and 
tan/? = 10): 



m d {Mqjjt) 



■ 0.94 MeV , 
0.346 MeV , 



m s {Mqut) 
mfj, {M GU t) 



■■ 17 MeV , 
73.0 MeV 



m b {Mqut) 
m T (Mqut) 



0.98 GeV 
1.25 GeV 



(Ai 



However, we did not include in our analysis the low-energy supersymmetric threshold corrections to 
the bottom quark mass, which can substantially modify the value of tti^Mgut) UM- These corrections 
take the form: 

m b = (1 + e b tan^) y b v d , (A.9) 

with the coefficient e b given by: 



£6 



2a 3 /iM 3 f 2 2 



3tt 



hi' 



f(M; 



3 I 



m r 



b R 



where / is a loop function defined by: 



/ '(ml, ml, mj) 



rln 



in; 



+ 



vt Mt 



167T 2 mr 



2 2 2 \ 

m,- ,m,- , 



TYIq 



m; 



mz 



mi 



rln 



mi 



mi 



ni' 



in z 



in , 



in; 



mi 



(A.10) 



(A.ll) 



The function / is of order one for superpartner masses between 100 GeV and 1 TeV, yielding typical 
values of e& ~ 2%. Due to the tan/3 enhancement, the threshold corrections can thus reach the 20% 
level for tan/3 = 10. This is enough to reach m b (MQUT) — 1-17 GeV, making it possible for the 
antisymmetric corrections to the GUT-scale mass relation M d = M e to account for the measured down 
quark and charged lepton masses. 

Indeed, very good fits of the charged lepton and down quark masses can be obtained once super- 
symmetric threshold corrections are taken into account. Since the quality of the fit does not depend on 
the precise values of the Xj, we restrict ourselves to the case of a single 45 vev in the (1, 1, 15) direction, 
i.e. we set x\ = xi = 23 = 3 in Eq. (A. 3). For each set of parameters (/Xi,Si) providing a good fit, Md 



and M e are determined as well as the mismatch matrix U m introduced in Section 4.2 Since only U m 
is needed for the computation of the final baryon asymmetry, we concentrate on this matrix from now 
on. Let us introduce the following parametrization: 



U„ 



?<PT 









v(ef 2 ,e^,e^,5 m ) 



3 #S 







(A.12) 



where V is a CKM-like matrix with three real angles and a complex phase. As Eq. (A. 3) contains more 



parameters than needed to fit the down-type fermion masses, we fix some of them, such as |ei| and 
the phase of £2. For |£i| <C (£2!, which roughly corresponds to |ei| < 0.001 GeV, good fits generally 
correspond to values of the 6™ 's of the order of the Cabibbo angle or smaller, for instance: 



9 m 
12 



0.3 



•jm 
'13 



0.1 



QUI 

'23 



0.35 . 



(A.13) 



For larger values of |£i| (|£i| > 0.001 GeV), one can obtain larger (1,2) and (1,3) mixing angles, e.g.: 

1 , fl^ ss 0.2 , 6% w 0.2 . (A.14) 



'12 



'23 



In our numerical study of leptogenesis, we use different choices for U m , combined with a non- vanishing 
high-energy or Majorana phase that we fix at ir/4. The resulting four sets of parameters are displayed 
in the table below: 
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set 


nm 
°12 


nm 
#13 


6*23 


5 m 




<t>? 




C 


4>t 


vr/4 


1 


1.07 


0.22 


0.21 


5.80 


3.21 


4.37 


5.86 


0.87 


6.16 




2 


1.07 


0.22 


0.21 


5.80 


3.21 


4.37 


5.86 


0.87 


6.16 


$2 


3 


0.28 


0.089 


0.37 


0.062 


3.15 


3.12 


6.03 


2.94 


6.19 




4 


0.17 


0.066 


0.29 


0.23 


3.14 


0.54 


0.015 


6.27 


0.0032 





To illustrate the influence of U m on the leptogenesis parameters, we plot in Fig. 11 the CP asymmetry 
and the washout parameter mi« for different values of O 7 ^. Part of the effect which can be seen is 

s on the right-handed neutrino masses. 



£1t 

due to the influence of the 9Yi ' 



Case +-+ 



Case H — + 



10" 



10" 



10" 




10 



12 



10 



13 



.14 



10 1 
V R (GeV) 



10 



15 



10 



16 
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> 
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10 1 
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Figure 11: E\ r and m lA1 as a function of vr in the (+,—,+) solution, for <f>2 — 7r /4 and different values of 
012 € [0,7r/4]. All other parameters in f/ m and high-energy phases are set to zero. The reference case #"2 = 
(Um — 1) is plotted in black. The other input parameters are as in Fig. [2} 
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